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We propose a new description of the SU (N) Yang-Mills theory on a lattice, which enables one to 
explain quark confinement based on the dual superconductivity picture in a gauge independent way. 
' This is because we can define gauge-invariant magnetic monopoles which are inherent in the Wilson 

loop operator. For SU (3) there are two options: the minimal option with a single type of non-AbeUan 
magnetic monopole characterized by the maximal stability subgroup H ~ U{2) — SU{2) x f/(l), and 
the maximal one with two types of Abelian magnetic monopoles characterized by the maximal torus 
subgroup H — U{1) X U{1). The maximal option corresponds to a gauge independent reformulation 
of the Abelian projection represented by the conventional MAG. In the minimal option, we have suc- 
cessfully performed the numerical simulation of the SU{3) Yang-Mills theory on a lattice. We give 
preliminary numerical results showing the dominance of the non-Abelian magnetic monopole in the 
string tension obtained from the Wilson loop in the fundamental representation, and the infrared dom- 
inance of a decomposed field variable for correlation functions after demonstrating the preservation 
of color symmetry which was explicitly broken by the conventional MAG. 
keywords: dual superconductivity, quark confinement, monopole dominance 
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1. Introduction 

The dual superconductivity which is believed as the promising mechanism for quark confinement 
is conjectured to occur due to the condensation of magnetic monopoles, just as the ordinary super- 
conductivity is caused by the condensation of the Cooper pairs [|l|][Q]. In the dual superconductor, the 
dual Meissner effect squeezes the color electric flux between a quark and an antiquark into a tube like 
region to form the hadronic string. The relevant data supporting the validity of this picture have been 
accumulated by numerical simulations especially since 1990 and some of the theoretical predications 
have been confirmed by these investigations; the infrared Abelian dominance, magnetic monopole 
dominance, center vortex dominance and non-vanishing off-diagonal gluon mass, which are the most 
characteristic features for dual superconductivity!^] [Q][|5|]. However, they are confirmed only in spe- 
cific gauges such as the maximal Abelian (MA), Maximal Center (MC) and Laplacian Abehan gauge, 
which break color symmetry. 

We give a new description of the Yang-Mills fields theory on a lattice, which is expected to give 
an efficient framework to explain quark confinement based on the dual superconductivity picture. The 
description enables us to extract in a gauge-independent manner the dominant degrees of freedom that 
are relevant to quark confinement in the Wilson criterion in such a way that they reproduce almost 
all the string tension in the Unear inter-quark potential. We have already given a new framework for 
the lattice SU (2) Yang-Mills theory as a lattice version of the Cho-Faddeev-Niemi-Shabanov (CFNS) 
decomposition in a continuum theory, and presented numerical evidences for its validity by performing 
numerical simulations. [|6|][^] 

The issue of generating magnetic monopoles in SU (N) Yang-Mills theory has been investigated 
so far under the MAG which breaks the original gauge group SU (N) into the maximal torus subgroup 
H = U (1)^^^ Then, MAG yields N — I types of Abelian magnetic monopole, in agreement with the 
observation due to the homotopy group: 7i2{SU{N)/U{lf-^) = 7ii{U (l)^-'^) = Z^-^ . Therefore, 
it tends to assume that magnetic monopoles of N — I types are necessary to cause the dual Meissner 
effect for realizing quark confinement. However, it is not yet confirmed whether or not N — I types 
of magnetic monopole are necessary to achieve confinement in SU {N) Yang-Mills theory. Rather, we 
have a conjecture that a single type of magnetic monopole is sufficient to achieve quark confinement 
even in SU (N) Yang-Mills theory, once it is defined in a gauge-invariant way. In fact, this scenario was 
originally proposed in based on a non-Abelian Stokes theorem for the Wilson loop operator[n]. 



For G = SU (3) there are two options: the minimal option with the maximal stability subgroup H = 
U{2) = SU{2) X ?7(1) is a new one (overlooked so far) on which we focus in this Talk, while the 
maximal one with H = U{\) x U{\) being equal to the maximal torus subgroup corresponds to a 
gauge-independent reformulation of the Abelian projection represented by the conventional MAG as 
reported in the lattice 2007 conference [P]. 



2. New variables on a lattice 

First, we summarize the result of a new description of the lattice SU {N) Yang-Mills theory given in 
Ref.|^ as an extension of SU (2) case. We wish to construct a lattice formulation in which an ordinary 
link variable Ux^^i G G = SU (N) is decomposed in a gauge-independent manner into two variables 
and Vx^ji, i.e., Ux,^ =Xx^^Vx,^, so that only the variable Vx,^ carries the dominant contribution for quark 
confinement in agreement with the dual superconductivity picture. For this purpose, new variables in 



9 



A new description of lattice Yang-Mils theory and non-Abelian monopoles as the quark confine rAkihko Shibata 



the relevant description are supposed to be transformed by a group element Q.^ G G as 

Ux,n ^ = a,u,^^al^^ , (2. la) 

Vx,M^VvV=n,V,,^al+^, (2.1b) 

= -^v,/J ~ ^-V^-V.M^l' (2-lc) 

where Vx^^ is defined as a link variable and transforms just like the original Yang-Mills link variable 
Ux,ii , while Xji^ is defined like a site variable representing a matter field and transforms according to the 
adjoint representation. It is important to see that decomposed variables are required to be transformed 
by the full SU (3) gauge group. In the conventional Abelian projection, on the contrary, the Vx,^ field is 
identified with an Abelian part and is supposed to be transformed as an Abelian field. 



The relationship between lattice variables and gauge variables in continuum theory [jlOp are given 

by 

Ux,ii = ^x^^iVx.ti = exp {-ig J dx^'A^ix)^ , (2.2a) 

Vx,fi = exp (^-ig J dx^\^{x)^ = exp {-ige\ ^{x + epL/l)) , (2.2b) 

^x,pL = exp (-/g£X^(x)). (2.2c) 

In order to obtain the new variables by the decomposition respecting the gauge transformation 
property given in the above, we consider the extended Yang-Mills theory, called the master Yang-Mills 
theory (see Figure |I|), by introducing a single type of the color field: 

h, :=0.vdiag(l/^,l/y3,-2/\/3)0l G5?7(3)/i7(2) =G/^, (2.3) 

such that it transforms according to the adjoint representation under an independent gauge transforma- 
tion ©v G G = SU (3). The decomposition U = VX is determined by solving the defining equation; 

Dl [y]h, := i {Vx,^K+n - ^xVx,n) = 0, (2.4a) 
tr(X,,^h,) = 0, (2.4b) 



where eq.( 2.4a ) represents that h(jc) is covariantly constant in the background Vx.^ and eq.( 2.4b ) means 
that Xf,^ has the vanishing H-commutative part. ' The solution of the defining equations are obtained 
uniquely for given Ux,^ and hx by way of a newly defined variable Vx.^i which does not belong to SU (3): 

2\/3 24 

yx,ll ■= Ux,n + {^xUx,ii + Ux,^iK+ll) + -J^xUx,^\ix+ll- (2.5) 

In fact, the SU (3) variable Vx.^ is obtained using the polar decomposition, together with another SU (3) 
variable Xv,^ as 

1 /3 

yv,M=}^(det\^)" , Vx^:={Vx,^Vl^y"^Vx,^, X,-,„ = (2.6) 

In order to obtain the equipollent theory (written in terms of the new variables) with the original Yang- 
Mills theory, we impose the reduction condition which plays the role of reducing the extended gauge 



'it is shown [ p^ using a non-Abelian Stokes tiieorem that a set of defining equations is obtained as a necessary and 
sufficient condition for the Wilson loop operator to be dominated by the decomposed variable Vx.^i in the sense Wc[f/v,^] — 
(const.)Wc[V,,p]. 
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Figure 1: The master Yang-Mills theory: 



symmetry SU (3)^ x [SU (3) /U {2)]q to the original gauge symmetry SU (3)^=0 so that the new vari- 
ables Vx,^,Xx^^ and hj^ transform under the same group Q. = & £ SU (3). Such a reduction condition 
can be given by minimizing the functional Frc of U_x,^ and h;^ under the independent local gauge trans- 
formations {n,0}: 

Frc[^,&;U.,h,K] = £tr ((D^p[/,,^]%)(D^p[/,,^]%)1-) . (2.7) 



By definition, the reduction condition (2.7) is invariant under the gauge transformation = and 
does not break the original gauge symmetry as expected. Therefore, we can impose any gauge fixing 
afterwards. 



3. Numerical simulations 

Next, we proceed to the numerical simulations. The link variable configurations {Ux^^} can be 
generated by using the standard method, since the decomposition Ux,^ = Vx^^Xx^^ is done in a gauge- 
invariant manner. The color field hj^ is determined so as to fulfill the reduction condition using the 
same algorithm as that used in the the gauge fixing procedure. We have generated the link variable 
configurations for the Wilson lattice action using 16'^ lattice at j8 = 5.70, and obtained the new variables 
by imposing the lattice Landau gauge (LLG) for the original Yang-Mills field, rather than the MA gauge 
which breaks color symmetry explicitly. 

We first focus on the color field. The new variables are defined so as to preserve the local gauge 
symmetry of the original Yang-Mills theory. Even after imposing LLG, therefore, the Yang-Mills 
theory should respect the global gauge symmetry, i.e., the color symmetry. This issue can be tested 
by measuring the vacuum expectation value (VEV) of the color field and its correlation functions. 
The VEV of a color field shows (/j'*) = for A = 1,2, ..,N^ — 1 = 8. In Figure ^, we give a plot of 
correlation functions (/j'^(jc)/j^(0)) for A,B = 1,...,8. All the diagonal parts (/j'^(x)/j'^(0)) have the 
same non- vanishing correlations, while all the off-diagonal parts {h^{x)h^{Qif) (A ^B) are vanishing. 
Thus, the correlation functions is of the form 

{h^{x)h^{0)) = 5^^^D{x), (3.1) 

which implies that the color symmetry is preserved. 

A 
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Figure 2: The correlation functions of the color field. (Left panel) diagonal parts, (Right panel) off-diagonal 
parts. 



Next we investigate the static interquark potential through the new variables. The Wilson loop 



operator in the fundamental representation in the continuum theory is rewritten into 1 1 1 ] 



Wc[A]=tr 



Pexpig£A^{x)dx^ /tr(l)= ^ t/Ml^lzexp J5''^F^v[V] 

dl^[^]zCxp{igiK,Zz) + igiJ,N^)}, (3.2) 

where we have defined K := d*F, J := dF, := 5*0j:A~^ and Ni, := 50eA"^ with the four- 
dimensional Laplacian, A = d6 + dd. Here @y, is the vorticity tensor defined by 0^^ = Jj^dS^^ {X{o))d{x 
X(a)) on the surface £ : dL = C spanned by the Wilson loop C. Note that V^(x) field in eq(|3.2[) is 
equivalent to the (x) reproduced from the new variable V^.^ in the continuum limit. Therefore, the 
magnetic monopole contribution to the Wilson loop average on a lattice is given by 



{Wc[U]) ^ {Wc[V]) = {exp (^/£i^,,^S,,^ j j , (3.3a) 

^^,H=L^'\^-^')U^iaPYdaSi^i^ + ^), dc,Sip=Jp, (3.3b) 

with the lattice magnetic monopole current given by 

:= := -argtr , (3.4) 

where dx denotes the forward difference (lattice derivative) in the A direction: dxf{x) := /(x + eX) — 
f{x). It should be noticed that the magnetic monopole current Kx^^ is gauge invariant, as can be 
seen from the transformation law of the new variables. This monopole should be identified with a 
"non-Abelian" magnetic monopole, since the monopole current is defined from Vv,/i which involves 
hv G G/H = SU {3)/U (2). The left panel of Figure || shows the distribution of the magnetic charge 
of the relevant magnetic monopoles obtained from 200 configurations, which have integral quantized 
values: rix^^ = Kx^^/ln € Z. The right panel of Figure || show the static potential calculated from 



the monopole part according to eq( 3.3a ). The numerical data of the static potential Vm{R) extracted 
from (exp {iT,x,fif^x,n'^x,n)) is well fitted by a function Vm{R) = —a,„/R + GniR with the value Om = 



0.13301(36). In comparison with the full string tension, for example, the result in Ref.[13], Ofuii 
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Figure 3: (Left panel) the magnetic monopole charge distribution, (Right panel) the static interquark potential 
calculated only from the magnetic monopole part. 

0.3879(39), thus, we have shown the non-Abelian magnetic monopole dominance for the string tension 
in the SU{3) Yang-Mills theory: 

a„,/a/„//= 0.87 ±0.19. (3.5) 



Finally, we devote to the correlation functions. The 2-point correlation functions (propagators) of 
the new variables and the original Yang-Mills gauge potential are defined by 



Dooix-y) := (C^(x)C^(y)> for Of,{y') e {V,.,^,Xy,^, Ay, 



(3.6) 



where an operator O^ix) is defined by the linear type, e.g., Ay,^ := — ul^)/2ge. The left panel 
of Figure ^ shows the behavior of correlation functions, Daa, Dyv and Dxx vs. the distance I = \x — y\. 
The correlation function Dyv damps slowly and has almost the same damping as Daa, while the Dxx 
damps quickly. This result suggests that the V part of the gluon propagator is dominated in the infrared 
region and mass generation by X part. As the variable Xx^^ transforms as an adjoint matter, (see 



eq.(2. Ic)), the Yang-Mills theory can have a gauge-invariant mass term = ^x^^O^ji)- The Fourier 
transformation of the massive propagator behaves for large Mxr as 



Dxx{r) 



-Mxr 



(3.7) 



(271)4 P+M^ 2(271)3/2 r3/2 ' 

and hence the scaled propagator r'/^Dxxir) should be proportional to exp(— Mxr). The right panel of 
Figure ^ shows the logarithmic plot of the scaled propagators in LLC This suggests the mass genera- 
tion of gluons, although there exists the mixing between V and X in LLG and Mx cannot be identified 
with the mass of X straightforwardly. 



4. Summary and discussions 

We have given the new description of the lattice Yang-Mills theory to give a gauge independent 
decomposition of the link variable, Ux,^ = Xx.^Vx,^ , which approves the the gauge independent formula- 
tion of the dual superconductivity picture. We have performed the numerical simulation in the minimal 
case of the SU (3) lattice Yang-Mills theory, and have demonstrated the color symmetry restoration, 
the non-Abelian monopole dominance in a gauge invariant way and the infrared V dominance in LLG 
other than MAG. We have also shown based on a non-Abelian Stokes theorem that the string tension 
for the fundamental quark is explained by the non-Abelian magnetic monopole defined from V part in 
the minimal case. 



A new description of lattice Yang-Mils theory and non-Abelian monopoles as the quark confine rAkihko Shibata 



1 [ , , , , , , , J 

! <AA> I — • — I 

•- <VV> ' ' 

J <XX> M« 

, <xv> 

0.1 - ■ B • 

* } " 

} 1 ' 

0.001 - ' 

0.0001 I ' ' ' ' ' ' ' 

'■""o 2468 10 12 14 16 012345678 

distance (I) distance (I) 

Figure 4: The correlation functions of original gauge field and new variables. (Left panel) The logplot of 
correlation functions. (Right panel) The scaled plot. 

To establish the dual superconductivity picture by the new description of the lattice Yang-Mills 
theory, we need further study for the full V dominance, e.g., "electric" and "magnetic" from the V part 
reproduces the full string tension, the N-ality property, and so on. These subjects are under investigation 
and will be discussed in a separate paper. 
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